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We model a social network by a random graph whose nodes represent agents and links between 
two of them stand for a reciprocal interaction; each agent is also associated to a binary variable 
which represents a dichotomic opinion or attribute. We consider both the case of pair-wise (p — 2) 
and multiple (p > 2) interactions among agents and we study the behavior of the resulting system by 
means of the energy-entropy scheme, typical of statistical mechanics methods. We show, analytically 
and numerically, that the connectivity of the social network plays a non-trivial role: while for pair- 
wise interactions (p — 2) the connectivity weights linearly, when interactions involve contemporary 
a number of agents larger than two (p > 2), its weight gets more and more important. As a 
result, when p is large, a full consensus within the system, can be reached at relatively small critical 
couplings with respect to the p = 2 case usually analyzed, or, otherwise stated, relatively small 
coupling strengths among agents are sufficient to orient most of the system. 

I. INTRODUCTION 

In recent years there has been an increasing interest in exploring a, quantitative approach 
to social sciences by means of statistical physics methods S, 0, Hi, 12, 13 1. In particular, 
a great attention has been paid to the topological structure representing the network of inter- 
action between agents, which is known to deeply influence the global behavior. Social systems 
are effectively envisaged by graphs, whose nodes represent agents and links between them rep- 
resent a couple interaction, and if relations between multiple agents are present, the network 
is represented by an hypergraph, formed by n-simplex interactions Q. In this framework, the 
behavior of a large number of interacting units, i.e. a social group, can be investigated by 
means of the energy-entropy scheme (free-energy variational principle), typical of statistical 
mechanics methods. 

The description of both the network topology, i.e. who is interacting with whom, and 
of the coupling strength between agents, i.e. the magnitude of each link, can be attained 
via a cost function H. In particular, we choose a class of such functions which may mimic 
several different contexts among N agents: Each cost function displays the mean-field form 
Ei!<i,< <i Jii h ■■■i P < J ii x °»a x ••■ x a i P > where crj = ±1, and the tensor J tunes the interaction 
strength in the network. Restricting for simplicity to the cases p = 2,3, we deal with agents 
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interacting in couples and triplets (strictly speaking, for p = 2 we talk about networks and for 
p > 2 we talk about hypergraphs [7|]). 

As for the network topology, many structures have been proposed along recent years, and 
they are in general built starting from three main architectures: the random graph, the small- 
world and the scale- free network which reproduce some general features of the observed 
social systems; here we focus on the former as it allows a simpler mathematical approach. 

A social system with this simple formal description can be used to model decision making: 
Oi is the opinion of the agent i and each agent i tries to aline his opinion (in the case of 
imitative behavior, i.e. Jij > 0) to other agents' viewpoint, interacting one by one (p = 2), or 
in larger groups (i.e. p > 3). 

Another appealing application of this models concerns trading among agents: Suppose we 
represent a market society only with couple exchanges (p = 2). Then, there are just sellers 
and buyers and they interact only pairwise. In this case if the buyer i has money (cr^ = +1) 
and the seller j has the product (<jj = +1), or if the buyer has no money and the seller has no 
products (pi = <7j = —1), the two merge their will and the cost function reaches the minimum. 
Otherwise, if the seller has the product but the buyer has no money (or viceversa), their two 
states are different (piOj — —1) and the cost function is not minimized. In this scenario, 
the possibility that an agent is satisfied increases only linearly with the number of his/her 
acquaintances, namely with the degree of the relevant node. In fact , the higher the number of 
"neighbors" , the larger the possibility of trading. 

When switching to the case p = 3, other strategies (on the timescale by which the connec- 
tivity remains constant) are available: for example the buyer may not have the money, but he 
may have a valuable good which can be offered to a third agent, who takes it and, in change, 
gives to the seller the money, so that the buyer can obtain his target by using a barter-like 
approach. In this case the contribution of the third agent k can either avoid the two frus- 
trated configurations of the previous picture, by providing a factor <7k = —1, or it can leave 
the frustration unaffected if he does not agree (er^ = +1). Interestingly, we find that in this 
case (p = 3), the amount of acquaintances one is in touch with (strictly speaking, the degree 
of connectivity a) does not contribute linearly as for p = 2, but quadratically: this seems to 
suggest that if a society deals primarily with direct exchanges, no particular effort should be 
done to connect people, while, if barter-like approaches are allowed, then the more connected 
the society, the larger the satisfaction reached on average by each agent in his specific goal. 
Intuitively, the above scenario seems to match the contrast among the classical barter-like 
approach of villages, where, thanks to the small amount of citizens, their degree of reciprocal 
knowledge is quite high and the money-mediated one of citizens in big metropolis, where a 
real reciprocal knowledge is fewer. 

In this work we want to pave both the analytical and the numerical analysis of what settled 
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so far: to tackle this task, at first we build our cost function in the next section 2, then 
in section 3 we analyze it by equilibrium statistical mechanics of disordered systems and in 
section 4 we corroborate our findings by Monte Carlo numerical simulations. 

II. DEFINITION OF THE MODEL 

First of all, we define a suitable Hamiltonian (a cost function) acting on a random network 
with connectivity a made up of N agents er,; = ±1, i £ [1, N]. 

Introducing p families of i.i.d. random variables {i*}, uniformly distributed on 

the previous interval, the Hamiltonian is given by the following expression 

H N (a,^(a)) = -J ^2 a i iC7 i 2...a i v (1) 

where, reflecting the underlying Erdos-Renyi graph, k is a Poisson distributed random variable 
with mean value -f(a)N and J > is the interaction strength, supposed to be the same for each 
p-plet. The relation among the coordination number a and 7 is 7 ~ a v . this will be easily 
understood a few lines later by a normalization argument coupled with the high connectivity 
limit of this mean field model. 

The quenched expectation of the model is given by the composition of the Poissonian average 
with the one performed over the families 

E[.] = E P EM = £ 6 - x ± £ [•]. (2) 

Following a statistical mechanics (SM) approach, we know that the macroscopic behavior of 
the system as a function of the average degree a and of the interaction strength J, is described 
by the following free energy density 

A(a,J)= lim A N (a, J) = lim -J-E In V exp(- J ff A r(cr, 7(a); J)). (3) 

(j 

The normalization constant can be extracted performing the expectation value of the cost 
function: 

00 -jN( ]\r\k j }- N 

nm = -E e fc V >< E E^^-^ - -iJNm^ ( 4) 

by which it is easy to see that the model is well defined and, in particular, it is linearly 
extensive in the volume. Then, in the high connectivity limit, each agent interacts with all 
the others and, in the thermodynamic limit, the coordination number a — > 00 as N. Now, if 
p = 2 the amount of couples in the summation scales as N(N — l)/2 and 7 = 2a provides the 
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right scaling; if p = 3 the amount of triples scales as N(N — 1)(N — 2)/3! and 7 = 3!a 2 again 
recovers the right connectivity behavior. The result ca be generalized to every finite p < N . 
Finally, we introduce the fundamental quantities expressed by the multi-overlap 

1 N 

1 (!) (™) (K\ 

qi...n = Jj2-^ a i - a i ' ( 5 ) 

i=l 

with a particular attention to the magnetization m = q\ = (1/N) E<=i a i- This plays the role 
of order parameter, representing the average opinion in decision making and the average trade 
in market. 



III. ANALYTICAL RESULTS 



In this section we summarize the scheme developed in statistical mechanics: Our goal is 
finding an explicit expression for the minimized free energy, which describes the overall behav- 
ior of our agents. To this task we decompose this quantity via the next equation ((6]) (whose 
proof is known in SM 0, into two quantities which can be estimated in an easier way, 
namely a cavity function <5 and a connectivity shift d a A: 



ry d 

A(a,J) =ln2 -—A(a,J) + *(a,J), 

p — 1 da 

where the cavity function \I/(a, J) = limjv^oo \I/jv(7, J) is defined, at finite N, as 

. k 2-i 
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Hence, we now need to evaluate the cavity function and the connectivity shift (the a derivative 
of the free energy density) . Starting with the latter and using the following properties of the 
Poisson distribution 

E[kg(k)]=N 7 E[g(k + l)], ^E[g(k)] = NE[g(k + 1) - g(k)}, (7) 



we can write 

d ., „ Id 



T,J) 
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(p l)a*- 2 E[ln JV Ej±i a *>" a * ~ ln£ * 



52J=i CT,! ■■■ cr i p 



Now, considering the following relation and definition 

e Ja^...<j iP = co S h J + o-j!... o-jp sinh J, 6* = tanhJ, 



(8) 
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Figure 1: Magnetization (main figure) and its normalized fluctuations (inset) for systems of different 
sizes and different dilution as a function of J q p_1 . Left panel: p = 2, right panel p — 3. The collapse 
of all the curves provides a strong evidence for the scaling of the interaction strength, coupled to the 
connectivity. 



and expanding the logarithm, we obtain 

-j-A(a, J) = (v— l)a p ~ 2 In cosh J - (p - l)a p ~ 2 - — >-6 n (ql _). (9) 
da n 

n—l 

With the same procedure, posing /?' = 2aO(J), and using a little of algebra [si, [fit, it is possible 
to show that 

J) = 2a^ J In cosh J + ^(m 2 ^) - ^!(^(p-D) + q^) 

and we get the result: a polynomial form for the free energy density 

A(a, J) = In 2 + a^ 1 lncosh J + y (^'(m 2 ^ 1 )) - (m p )) + 0(^' 3 ). 

It is straightforward to see that for p = 2 the well known diluted Curie- Weiss is recovered as 
well as its criticality at (3' = 1 Further, it is enough to explore the ergodic phase (m = 0) 
to see that a appears at the power p — 1 (i.e. 2 for p = 3). 



IV. NUMERICAL RESULTS 



We now analyze the system introduced in Sec. 2, from the numerical point of view by 
performing extensive Monte Carlo simulations [3| ; here we especially focus on the cases p = 2 
and p = 3. In general, we find that the system is able to relax to a well-defined steady 
state characterized by average observables (such as total energy and magnetization) which are 
independent of the initial configuration and of the system size N, as long as N is large enough 
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to avoid finite-size effects. On the other hand, the average observables vary as the average 
coordination number a and/or the interaction strength J are tuned. In particular, as shown 
in Fig. 1 we find that the curves for the order parameter collapse when plotted as a function 
of J a p ~ l , confirming the scaling found analytically. We also notice that the system exhibits a 
phase transition at a critical interaction strength J c , which depends on the connectivity of the 
underlying network according to J c = or 1 for p = 2 and J c = Aa~ 2 for p = 3. The existence 
of a phase transition is also confirmed by the corresponding peak displayed by the fluctuations 
of the order parameter (to 2 ) — (to) 2 , see insets in Fig. 1. 

Before concluding, we stress the non-trivial role played by the connectivity of the social 
network: while for pair-wise interactions (p = 2) a weights linearly, when interactions involve 
contemporary a number of agents larger than two (p > 2), its weight gets more and more 
important. As a result, whenp is large, a full consensus within the system, i.e. a ferromagnetic 
state, can be reached at relatively small critical strengths with respect to the p = 2 case usually 
analyzed, or, otherwise stated, relatively small coupling strengths among agents are sufficient 
to orient most part on the social system Q . 
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